The percolation transition of hydration water: from planar hydrophilic surfaces to 

proteins 
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The formation of a spanning hydrogen-bonded network of hydration water is found to occur via 
a 2D percolation transition in various systems: smooth hydrophilic surfaces, the surface of a single 
protein molecule, protein powder and diluted peptide solution. The average number of water-water 
hydrogen bonds uh at the percolation threshold varies from 2.0 to 2.3, depending on temperature, 
system size and surface properties. Calculation of nn allows an easy estimation of the percolation 
threshold of hydration water in various systems, including biomolecules. 

PACS numbers: 61.20Ja, 64.60Ak, 87.15Aa 



The existence of a spanning network of hydration 
water in biosystems enables their biological functions 
IH S [3 S With increasing hydration level, an 
ensemble of finite (non-spanning) clusters of hydration 
water transforms via a quasi-2D percolation transition 
into a state with a spanning water network 0,0. The 
full dynamics of biomolecules is restored, when they are 
covered by about a "monolayer" of water. The first ap- 
pearance of such a "monolayer" corresponds to a quasi- 
2D percolation transition of the hydration water at the 
surface of a single biomolecule 0, HI ■ Simulation studies 
of various properties of hydrated biosystems below and 
above the percolation threshold of the hydration water 
can help to clarify the role of the spanning water net- 
work in the onset of biological functions. Such studies 
require the knowledge of the percolation threshold of wa- 
ter in the system of consideration. This information can 
be obtained by conventional percolation analysis of wa- 
ter clustering 0, 0, 0, which is extremely time 
consuming. In the present paper we propose a simple 
method to locate the percolation threshold of hydration 
water even in complex systems, using the average number 
of water- water hydrogen bonds. This method is derived 
from extensive computer simulation studies of the perco- 
lation transition of hydration water in various systems: 
water adsorbed at smooth hydrophilic planes and spheres 
at surfaces of rigid and flexible single lysozyme 
molecules and lysoz yme powder , hydration water 
in protein solutions 

Water molecules are considered to belong to the same 
cluster if they are connected by an uninterrupted path 
of hydrogen-bonds. Water clustering was studied by 
computer simulations at various hydration levels in low- 
hydrated systems 0, Q] and at various temperatures 
and widths of the hydration shell for proteins in solution 
[l5j |. The percolation threshold of water at planar sur- 
faces and in a model lysozyme powder (both systems are 
infinite due to the periodic boundary conditions) , can be 
located by using the spanning probability, cluster size dis- 
tribution and fractal dimension of the largest cluster. In 
closed systems, such as a surface of a finite object, an in- 



finite cluster can not appear and the spanning probabil- 
ity is not defined. However, the formation of a spanning 
network of hydration water at the surface of a smooth hy- 
drophilic sphere and of a single biomolecule, was found 
to occur in a way similar to infinite systems |7llS U13l , ll5j . 




FIG. 1: Cluster size distribution ns in the hydration shell 
of an elastin-like peptide at T = 300 K for several widths 
D of the hydration shell, which corresponds to the following 
numbers N w of water molecules in the shell: 135, 141, 147, 
153, 159 and 164 (from bottom to top). Power laws for 2D 
and 3D percolation thresholds are shown by solid and dashed 
lines, respectively. 

The probability distribution ns of cluster sizes S obeys 
a universal power law ns ~ S~ T in a widest range of 
S at the percolation threshold, with r rs 2.05 for 2D 
and r m 2.2 for 3D percolation HJ. In Fig. 1 we show 
the cluster size distributions ns of water molecules in 
the hydration shell of the smallest studied system, the 
fully hydrated elastin-like peptide GVG(VPGVG) 3 , for 
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various sizes of the shell width D. A hump in ns at large 
S reflects the truncation of the large clusters due to the 
finite size of the hydration shell. The behavior of rig 
allows the location of the percolation threshold between 
N w = 147 and N w = 153, even without any assumption 
about the dimensionality of the transition. 

The dimensionality of the percolation transition can be 
determined from the effective fractal dimension df of the 
largest water cluster, which can be obtained from a fit of 
its mass distribution: m(r) ~ r df . The range of r, used 
in the fits of m(r), is determined by the system size. We 
limited this range by the half of the simulation box L for 
planar surfaces and lysozyme powder, by the diameter of 
the spherical surfaces and by the shortest value of the ef- 
fective axis of a single lysozyme molecule, which is about 
24 A. At the percolation threshold, the fractal dimension 
of the largest cluster is dj D 1.9 for a 2D system and 
d^ D » 2.53 for a 3D system ||. We have found that in 
low-hvdrated systems the value of df is close to the value 
df D 0, EE Q- Tnis indicates the quasi-2D character 
of the percolation of adsorbed water. In particular, the 
threshold hydration level and the dimensionality of the 
water percolation transition in model lysozyme powder 
remarkably agrees with experimental results 4J. Note 
that the fractal dimension of the largest water cluster in 
the hydration shell of the elastin-like peptide in solution 
can not be determined due to the small system size (the 
average radius of gyration of this peptide does not exceed 
about 8 A). 

TABLE I: Number of water molecules in the first (Ni) and sec- 
ond (N2) hydration shell per unit surface area and the average 
number of water-water hydrogen-bonds nu at the quasi-2D 
percolation thresholds of various model systems and proteins. 



System 


T 1 K 


Ni / 100 A 2 


N 2 1 Ni 


nu 






(± 0.05) 


(± 0.005) 


(± 0.01) 


Plane 80x80 A 2 


425 


7.5 


0.043 


2.20 


Plane 100x100 A 2 


425 


7.4 


0.043 


2.22 


Plane 150x150 A 2 


425 


7.3 


0.044 


2.22 


Sphere R = 15 A 


425 


8.6 


0.12 


2.11 


Sphere R = 30 A 


425 


8.5 


0.11 


2.14 


Sphere R = 50 A 


425 


8.3 


0.10 


2.15 


Sphere R = 15 A 


475 


9.1 


0.36 


1.95 


Lysozyme powder 


300 


2.1° 




2.32 


Rigid lysozyme 


300 


4.9 


0.38 


2.30 


Flexible lysozyme 


300 


5.0 


0.44 


2.30 


Lysozyme powder 


400 


2.5° 




2.02 


Rigid lysozyme 


400 


5.1 


0.98 


2.05 


Elastin 


320 


10.8 




2.08 


Elastin 


300 


10.3 




2.10 


Elastin 


280 


9.8 




2.10 


Elastin 


260 


9.2 




2.08 



"total number of water molecules per one lysozyme 



The percolation thresholds of hydration water in all 



systems studied are given in Table I in terms of N\, the 
number of water molecules in the first hydration shell per 
unit surface area. For lysozyme powder, N\ is simply the 
total number of water molecules in the simulated system 
divided by the number of protein molecules. The values 
obtained for the threshold hydration level in terms of Ni 
noticeably vary from one studied system to another and 
also depend on temperature. These variations of Ni at 
the threshold should be expected, as this surface cover- 
age is similar to the occupancy variable in lattice models, 
which is essentially non-universal. For random percola- 
tion in 2D lattices, the threshold value of site and bond 
occupancy depends on the lattice type |IJ. The honey- 
comb lattice with 3 neighbors and the square lattice with 
4 neighbors seem to be the most relevant for water ad- 
sorbed at hydrophilic surfaces [lilllijl . For these lattices, 
the threshold values of the occupancy variable are ~ 0.70 
and ~ 0.59 for site percolation and ~ 0.65 and 0.50 for 
bond percolation, respectively. These threshold values 
are more universal in terms of the average number of 
bonds, which is ~ 2.09 and ~ 2.37 for site percolation 
and ~ 1.96 and 2.00 for bond percolation, on the hon- 
eycomb and square lattice, respectively. As the water 
percolation process can be considered as a correlated site- 
bond percolation [l^ , one may expect a rather universal 
threshold value in various hydrated systems in terms of 
Tin, the average number of water- water H-bonds, formed 
by each water molecule. Such expectation is supported 
by the observed similarity between percolation in bulk 
liquid water and random bond percolation in 3D lattice. 
In bulk liquid water, nu at the percolation threshold is 
about 1.53 i- e -i rather close to the value 1.55 for 
bond percolation in a 3D diamond lattice with 4 neigh- 
bors Q- 

The dependence of the fractal dimension df on nn for 
some of the systems studied is shown in Figs. 2 and 3. 
Below the percolation threshold, df is essentially an effec- 
tive fractal dimension, because most of the largest water 
clusters are not fractal objects. That is why the values 
of df noticeably depend on the system size and geometry 
at low hydration levels. At the percolation threshold, the 
structure of the largest water cluster is close to a fractal 
and df approaches the threshold fractal dimension dj D 
for 2D percolation. For three studied planar surfaces 
of different sizes, the percolation thresholds of hydration 
water practically coincide: at T = 425 K they are located 
at n H = 2.21±0.01 and N x = 0.074±0.001 A~ 2 . At the 
spherical surfaces, the percolation threshold corresponds 
to Ni = 0.085±0.002 A~ 2 , that is about 15% higher than 
the threshold value of N\ for planar surfaces. However, 
the value of nu at the threshold is 2.13±0.02 for spher- 
ical surfaces, that is just 4% lower than for the planar 
surface. The higher threshold hydration levels N\ for the 
spheres are accompanied by higher concentrations of wa- 
ter molecules in the second hydration shell N2 . The ratio 
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FIG. 2: Fractal dimension of the largest cluster dj as a func- 
tion of the average number uh of H-bonds between water 
molecules near smooth hydrophilic planar (lower panel) and 
spherical (upper panel) surfaces. 



N2/N1 is almost 12% for spheres, whereas it is only 4% 
for the planar surfaces of the same hydrophilicity (see Ta- 
ble I) . The presence of water molecules in the second shell 
effectively reduces the number of water- water H-bonds in 
the first shell, thus resulting in a higher threshold value 
of Ni and a lower threshold value of njj. Such shift of 
tih to lower values can be considered as a trend toward 
three-dimensionality, where tih is about 1.53 at the per- 
colation threshold of bulk liquid water . 

A similar trend toward larger threshold values of N\ 
and smaller values of is observed with increasing tem- 
perature. When the temperature rises from T = 425 to 
475 K, the threshold value of Ni for a sphere with radius 
R = 15 A increases from 0.086 to 0.09lA , whereas 
the fraction of molecules in the second hydration shell 
increases and N2/N1 changes from 0.12 to 0.36. This 
is accompanied by a decrease of uh from 2.11 to 1.95. 
Hence, an increase of temperature by 100 degrees causes 
a decrease of the nn value at the threshold by about 15 
%. This trend corresponds to the growing importance of 
the "bond" percolation relatively to the "site" percola- 
tion with increasing temperature in site-bond percolation 
of water. 

At the surfaces of the rigid lysozyme molecule, the 
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FIG. 3: Fractal dimension of the largest cluster Aj as a func- 
tion of the average number uh of H-bonds between water 
molecules at the surface of rigid (open squares) and flexible 
(solid squares) lysozyme and in the hydrated lysozyme pow- 
der (open circles). 



flexible lysozyme molecule, and in the model lysozyme 
powder, the percolation transition of hydration water at 
ambient conditions (T — 300 K) occurs when the aver- 
age number of water-water H-bonds uh is 2.31±0.01 (see 
Table I and Fig. 3). Intriguingly, the value of hh at the 
percolation threshold remains highly universal even in 
such complex systems as hydrated proteins. An increase 
of the temperature to T — 400 K reduces the threshold 
value of nn to 2.03±0.02, i.e by about 15 %. Hence, 
the effect of temperature on the threshold value of uh 
is remarkably similar at protein surfaces and at smooth 
surfaces of essentially hydrophilic spheres (compare the 
values N2/N1 for these systems in Table I). 

The surface of a biomolecule in dilute aqueous solution 
is completely covered by hydration water, which forms a 
spanning network at low temperature. Upon heating, 
this spanning water network breaks up via a percola- 
tion transition, and the hydration water shell becomes 
an ensemble of finite (non-spanning) water clusters |l5| . 
The location of the percolation threshold at the surface 
of fully hydrated elastin-like peptide, using cluster size 
distribution n$, indicates, that at T = 300 K, nu is 
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about 2.09 ± 0.02 at the threshold (see Fig. 1). More 
accurately, the percolation threshold can be located by 
additional analysis of t he p robability distribution of the 
largest cluster size [1J, [l5| . The estimated temperature 
of the percolation transitions of hydration water in the 
shells of various width D are shown by open circles in 
Fig. 4 as a function of the average value of water-water 
H-bonds uh- Evidently, the value of njj at the percola- 
tion threshold is about 2.1 for any reasonable choice of 
the hydration shell width D. The value nu for hydration 
water in protein solutions decreases almost linearly with 
temperature (see Fig. 4 and also Ref. This gives us 

the possibility to estimate the value of nu in a wide tem- 
perature range by simulations at two temperatures and, 
subsequently, to estimate the temperature of the thermal 
breaking of the spanning network of hydration water. 
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FIG. 4: The average value nn of water-water H-bonds in 
the hydration shell of an elastin-like peptide as a function of 
temperature (closed circles and dashed lines). The widths of 
the hydration shell D are: 4.40, 4.45, 4.50, 4.55, 4.60, 4.70, 
4.80, 4.90, and 5.00 A. Percolation thresholds are shown by 
open circles. 

Summarizing our studies of the quasi-2D percolation 
transition of hydration water in various model systems, 
we propose a relatively simple method to indicate the ex- 
istence of a spanning network of hydration water by the 
analysis of the average number nu of hydrogen bonds 
between the water molecules in the hydration shell. This 
value can be obtained by conventional computer simula- 
tions or estimated from experimental data |2fJ. At am- 
bient temperature, the threshold value of uh is about 2.1 
for fully hydrated systems and about 2.3 for low-hydrated 



systems. A value uh above the threshold indicates the 
presence of a spanning network of hydration water. The 
threshold value of n H decreases slightly with increasing 
temperature. 
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